n-p Interaction Effects on the Double Beta Decay Nuclear 
Matrix Elements for Medium Mass Nuclei 



P. K. Raina and A. Shukla 
Department of Physics and Meteorology, IIT Kharagpur-721302, India 

P. K. Rath, B. M. Dixit, K. Chaturvedi, and R. Chandra. 
Department of Physics, University of Lucknow, Lucknow-226007 , India 

S. K. Dhiman and A. J. Singh. 

Department of Physics, H. P. University, Shimla-171002, India. 

Abstract 

The quality of HFB wave functions are tested by comparing the theoretically calculated 
results with the available experimental data for a number of spectroscopic properties like 
yrast spectra, reduced B(E2) transition probabilities, quadrupole moments and g-factors for 
the nuclei involved in 2v (3(3 decay. It is observed that the np interactions vis-a-vis the 
deformations of the intrinsic ground states of medium mass nuclei play a crucial role in the 
fine tuning of the nuclear matrix elements, M21, . 

It is well established by now that the implications of nuclear (3(3 decays are far reaching in 
nature in general. Ov (3(3 decay in particular is one of the very rare promising processes to test the 
physics beyond the standard model (SM) of fundamental particles. These aspects of nuclear (3(3 
decay have been excellently elaborated in a number of review articles over the past years [1-8]. 

The 2v (3(3 decay, a second order process of weak interaction that conserves the lepton number 
exactly, is allowed in the SM. The half life of 2v (3(3 decay is a product of accurately known phase 
space factor and appropriate nuclear transition matrix element M2,,. The half lives of 1v (3(3 decay 
have been already measured for about ten nuclei and the values of can be extracted directly. 
Consequently, the validity of different models employed for nuclear structure calculations can be 
tested by calculating the M2,,. 

It is observed that in all cases the 2v (3(3 decay matrix elements are sufficiently quenched. The 
main motive of all the theoretical calculations is to understand the physical mechanism responsible 
for the suppression of the M-2v The M.2 V is calculated mainly in three types of models. One is 
the shell model and its variants. The second is the quasiparticle random phase approximation 
(QRPA) and extensions their of. The third type of models is classified as the alternative models. 
The details about these models - their advantages as well as shortcomings - have been discussed 
excellently by Suhonen and Civitarese [5] and Faessler and Simkovic [6]. 

All the nuclei undergoing (3(3 decay are even-even type. Hence the pairing degrees of freedom 
play an important role. Moreover, it has been conjectured that the deformation can play a crucial 
role in (3(3 decay rates. Hence it is desirable to have a model that incorporates the pairing and 
deformation degrees of freedom on equal footing in its formalism. For this purpose the Projected 
Hartree-Fock-Bogoliubov (PHFB) model is one of the most natural choices. Coincidentally, most 
of the (3(3 decaying nuclei fall in medium mass region. The success of the PHFB model in explaining 
the observed experimental trends in this mass region has motivated us to apply the PHFB wave 
functions to the study of nuclear (3(3 decay as well. 

The mass region A«100 provides us with a nice example of shape transitions [9], where, at one 
end, nuclei can be described in terms of shell model wave functions involving a small number of 



configurations and, at the other end of this region, we find good evidence of rotational collectivity. 
These nuclei lie between doubly magic 132 Sn and the strongly deformed 100 Zr, near which the 
structural changes are rather rapid with the addition of protons and neutrons. In the past there 
have been many attempts, [10-15] to explore the factors responsible for the structural changes in 
this mass region. 

Federman and Pittel [14] computed the deformation energy in the framework of Hartree-Fock- 
Bogoliubov (HFB) theory in conjunction with the surface delta interaction (SDI), suggesting that 
the neutron-proton (n-p) interaction in the spin-orbit partner (SOP) orbits-lg 9 /2 and lg7/2 in 
this case may be instrumental vis-a-vis the onset of deformation in Mo isotopes with A>100. 
A systematic study of the behavior of the low-lying collective states of neutron-rich even Cd, 
Pd, Ru, and Mo isotopes has lead to the conclusion that these structural changes are related to 
the exceptionally strong n-p interaction in this region. It has also been observed that the n-p 
interactions among the SOP orbits have a deformation producing tendency and the systematics of 
low-lying states are intricately linked with the nature of n-p interaction. 

The sensitivity of the yrast spectra and the transition charge densities (TCD) to the neutron- 
proton interaction strength has lead to the fixing of these strengths very accurately and has been 
demonstrated [15] through the examples of 110 Cd and 114 Cd. We have adopted this method for 
fixing the p-n strength of QQ interaction by looking at the spectra of 2 + state of the nuclei involved 
in (3(3 decay. 

A large number of theoretical as well as experimental studies of 2v (3(3 decay have already been 
done for /T/T of 96 Zr, 100 Mo, 110 Pd, 12 8,i30 Te nuclci and e + pp ( p+^Q and E CEC) in 

case of 96 Ru, 106 Cd, 124 Xe and 130 Ba nuclei over the past few years with more emphasis on 100 Mo 
and 106 Cd cases. The (3(3 decay is not an isolated nuclear process. The availability of data permits 
a rigorous and detailed critique of the ingredients of the microscopic models used to provide a 
description of these nuclei. 

We have studied the 2v (3(3 decay not isolatedly but together with other observed nuclear phe- 
nomena. This is in accordance with the basic philosophy of nuclear many body theory, which is to 
explain all the observed properties of nuclei in a coherent manner. Hence as a test of the reliability 
of the wave functions, we have calculated the yrast spectra, reduced B(E2) transition probabilities, 
static quadrupole moments and g-factors and compared with the available experimental data. 

The theoretical formalism to calculate the half life of 2v (3(3 decay mode has been given by 
Haxton and Stephenson [1] , Doi et al [2,3] and Tomoda [4]. Very brief outlines of the calculation 
of nuclear transition matrix elements of the (3(3 decay in the PHFB model are presented here. 
Details of expressions used in calculation of spectroscopic properties can be found in Dixit et al 
[16]. 

The half-life of 2v (3(3 decay for + -^0 + transition is given by 
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T?? 2 (0+^0+) =G 2 u\M 2 A (1) 



where 



and the integrated kinematical factor Q 2v can be calculated with good accuracy [8]. If the Ejy of 
Eq.(ll) is replaced by an average (En), the summation over intermediate states can be completed 
using the closure approximation one obtains 

2M% T 2M& 
(E N ) - (Mj + M F ) /2 E d W 

where the double Gamow- Teller matrix element (DGT) M"q T is defined as follows 
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Employing the HFB wave functions, one obtains the following expression for the (3(3 decay 
nuclear transition matrix clement. 
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and 



n { z,N U z-2,N + 2)(6) = {det[l + Fg^/W^]} 1 ^ x {det[l + ^^W]} 172 (7) 

The tt (v) represents the proton (neutron). PHFB calculations are summarized by the coeffi- 
cients (Vi m , Vi m ) an d Cij i7n and their matrices (Fjv ! z(0)) Q/ 3and ({n,z)- The details can be found 
in reference [16]. 

In the present calculations we treat the doubly even nucleus 76 Sr (Z=N=38) as an inert core 
and the valence space is spanned by the orbits lpi/2, 2si/2, ld 3 / 2) lds/2; 0g7/ 2 > 0§9/2 an( i Ohn/2 
for protons and neutrons. The set of single particle energies (SPE's) but for the e(0h n / 2 ) which 
is slightly lowered, employed here is same as used in a number of successful shell model as well 
as variational model [10-16] calculations for nuclear properties in the mass region A=100. The 
effective two-body interaction is the PPQQ type [17]. 

Table 1. Variation in excitation energies in MeV of J 7r =2 + , 4 + an d 6 + yrast states for 
100 Mo and 100 Ru nuclei with change in x pn keeping fixed G p = -0.30 MeV, G n = -0.20 MeV and 
e(0h n/2 ) = 8.6 MeV. 



Nucleus 


Xpn 


0.01826 


0.01866 


0.01906 


0.01946 


0.01986 


Exp. [21] 




E2+ 


0.6865 


0.5851 


0.5356 


0.4493 


0.3923 


0.5355 


100 Mo 


E4+ 


1.7028 


1.5333 


1.4719 


1.3070 


1.1861 


1.1359 




E6+ 


2.9355 


2.7213 


2.6738 


2.4560 


2.2854 






Xpn 


0.01758 


0.01798 


0.01838 


0.01878 


0.01918 






E2+ 


0.6597 


0.5923 


0.5395 


0.4930 


0.4445 


0.5396 


100 Ru 


E4+ 


1.8175 


1.6733 


1.5591 


1.4531 


1.3372 


1.2265 




E6+ 


3.2746 


3.0615 


2.8940 


2.7329 


2.5519 


2.0777 



The strength of the pairing interaction is fixed through the relation G p = -30/A MeV and G n = 
-20/ A MeV. These values of G p and G n have been used by Heestand et al [18] to successfully 
explain the experimental g(2 + ) data of some even-even Ge, Se, Mo, Ru, Pd, Cd and Te isotopes in 
Greiner's collective model [19]. The strengths of the like particle components of the QQ interaction 
are taken as: x PP = Xnn = -0-0105 MeV b~ 4 . These values for the strength of the interaction are 
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comparable to those suggested by Arima on the basis of an empirical analysis of the effective 
two-body interactions [20]. 

Table 2. Experimental half-lives T^j 2 , and corresponding nuclear matrix elements M 2v , along 
with the theoretical values in different models for + — > + 2vf3f3 decay of 100 Mo. The numbers 
corresponding to (a) and (b) are calculated for gA=l-25 and 1.0 respectively. 



Experiment 


Theory 


Ref 


Projects 


-4/2 

(10 18 yrs) 


\M 2v \ 


Ref 


Models 


\M 2v \ 


1/2 

(10 18 yrs) 
a) b) 


[22] 


UC-Irvin 


6.82±%i% ± 0.68 


aju.izo_ 009 


[16] 


PHFB 


0.152 


4.57 11.15 








b)0.195l» : ^ 


[24] 


SRPA(WS) 


0.059 


30.4 74.3 


[23] 


NEMO 


9.5±0.4±0.9 


a)0.106lg;88? 


[25] 


SU3(SPH) 


0.152 


4.59 11.2 








b)o.i65±g:8ig 


[25] 


SU3(DEF) 


0.108 


9.09 22.2 


[8] 


Average 


8.0±0.7 




[26] 


QRPA(EMP) 


0.197 


2.73 6.67 



As an illustrative case we look into the details of calculations for double beta decay of 100 Mo 
nucleus. The x pn ls varied so as to obtain the spectra of 100 Mo and 100 Ru in optimum agreement 
with the experimental results. In Table 1, we have presented the theoretically calculated yrast 
energies for levels of 100 Mo and 100 Ru for different values of x pn - ^ is clearly observed that as the 
X pn is varied by 0.0016 MeV b~ 4 , the E 2 , decreases by 0.2942 MeV in case of 100 Mo and 0.2152 
MeV in case of 100 Ru respectively. This is understandable as there is an enhancement in the 
collectivity of the intrinsic state with the increase of |x pl J, hence the E 2 decreases. The optimum 
values of Xpn corresponding to 100 Mo and 100 Ru are 0.01906 MeV b" 4 and -0.01838 MeV b~ 4 
respectively. Thus for a given model space, SPE's, G p , G„ and x PP -> we have fixed x pn through the 
experimentally available energy spectra. 

Table 3. Experimental limit on half-lives (T^j 2 ) and corresponding extracted matrix elements 

M 2v along with their theoretically calculated values for 2v [3 + [3 + \[3 + EC\ECEC decay of 106 Cd 
for + -^0 + transition. The numbers corresponding to (a) and (b) are calculated for gA=l-25 and 
1.0 respectively. 



Decay 
Mode 1, 


Experiment 


Theory 


Ref Tf /2 (yrs) 


Ref Model \M 2v \ (Tf /2 Y 
a) b) 


[3 + !3 + 
10 25 yrs 


[27] >2.4xl0 20 ** 
[28] >1.0xl0 19 * 
[29] >9.2xl0 17 


[31] PHFB 0.238 35.42 89.56 
[32] QRPA 0.166 72.79 # 180 # 
[28] QRPA 1.226 1.33 3.3 


f3 + EC 
10 21 yre 


[27] >4.1xl0 2 ° 
[28] >0.66xl0 19 * 
[29] >2.6xl0 17 


[31] PHFB 0.238 8.97 22.69 
[32] QRPA 0.169 17.79 # 44.0 # 
[33] SU(4) 0.198 13.00 32.15 


ECEC 
10 20 yrs 


[30] >5.8xl0 1Y 


[31] PHFB 0.238 11.24 28.42 
[32] QRPA 0.169 22.24 # 55. 0* 
[33] SU(4) 0.193 17.00 42.04 



* denotes half-life limit for Ov + 2v mode, ** denotes half-life limit for Ov + 2v +0vM mode, 
* shows half-life with WS potential, ^Number below the mode is multiplication factor of Half- life for 
theoretical values. 

From the overall agreement [16] between the calculated and observed electromagnetic prop- 
erties, it is clear that the PHFB wave functions of 100 Mo and 100 Ru generated by fixing x pn to 
reproduce the yrast spectra are quite reliable. 

The double beta decay of 100 Mo ^ 100 Ru for + — > + transition has been investigated by many 
experimental groups [22,23] as well as theoreticians by employing different theoretical frameworks 
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[24-26]. In Table 2, we have compiled some of the latest available experimental and the theoretical 
results along with our calculated M 2 „ and the corresponding half-life T^y 2 . We have used a phase 
space factor G 2 „ = 9.434 xlO~ 18 yr _1 given by Doi et al [2] and an energy denominator = 
11.2 MeV given by Haxton et al [1]. In column 4 of Table 2, we have presented the M 2 „ extracted 
from the experimentally observed T^ 2 using the phase space factor given above. The phase space 
integral has been evaluated for g y j=1.25 by Doi et al [2]. However in heavy nuclei it is more 
justified to use the nuclear matter value of gA around 1.0. Hence, the experimental M 2 „ as well as 
the theoretical Ty 2 are calculated for g J 4=1.0 and 1.25. The present calculation and that of Hirsch 
et al using SU3(SPH) [25] give nearly identical value. They are close to the experimental result 
given by De Silva et al [22] for gA=l- 25 while for g^=1.0, the above two M 2 „ are in agreement 
with the results of NEMO. The calculated values given by Stoica using SRPA(WS) [24] are too 
low and those from Suhonen et al [26] are slightly on higher side. Further the value M 2 „ given by 
Hirsch ct al using SU3(DEF) [25] favors the results of NEMO [23] for g^=1.25. 

Another example we take for the + -^0 + positron (3(3 {(3 + (3 + , /3 + EC and ECEC) decay of 
106 Cd^ 106 Pd. This transition has also been investigated by many experimental groups and in 
different theoretical frameworks. In Table 3, we have compiled some of the latest available experi- 
mental [27-30] and the theoretical results [31-33] along with our calculated M 2j/ and corresponding 
half-lives T^ 2 . We have used phase space factors given by Doi and Kotani [3] and the average 
energy from Haxton and Stephenson [1] . Our calculated values are nearly half of the recently given 
QRPA results of Suhonen and Civitaresse [32] for all the three modes. The theoretical values of 
PHFB and SU(4) [33] are in better agreement (factor of roughly two third) for the /? + EC and 
ECEC modes. 

From the above discussions, it is clear that the validity of nuclear models presently employed to 
calculate the M 2 „ cannot be uniquely established due to error bars in experimental results as well 
as uncertainty in g^. Further work is necessary both on the experimental and theoretical front to 
judge the relative applicability, success and failure of various models used so far for the study of 
double beta decay processes. 

As an example to see quantitatively the effect of deformation on M 2 „ vis-a-vis the variation 
of the strength of pn part of the QQ interaction, the results are displayed in Fig.l for 100 Mo and 
106 Cd cases. It is observed that the deformations of the HFB intrinsic states play an important 
role in the calculations of M 2 „ and hence on the half life. 

To summarize, we have first tested the quality of HFB wave functions by comparing the the- 
oretically calculated results for a number of spectroscopic properties of nuclei involved in double 
beta decay. To be more specific we have computed the yrast spectra, reduced B(E2) transition 
probabilities, quadrupole moments and g-factors. Some of the results have been presented for two 
very widely studied cases of (3~(3' decaying 100 Mo and e + (3(3 {(3 + (3 + , /3+EC and ECEC) decaying 
106 Cd nuclei. Reliability of the intrinsic wave functions for calculation of 2v (3(3 nuclear matrix 
elements, M 2 „ , has been discussed. Further, we have shown that the np interactions viz a viz 
the deformations of the intrinsic ground states of 100 Mo 100 Ru, 106 Cd and 106 Pd play important 
role in arriving at the appropriate nuclear matrix elements. A reasonable agreement between the 
calculated and observed spectroscopic properties as well as the 2v (3(3 decay rate of most of the 
nuclei in medium mass region makes us confident to employ the same PHFB wave functions for 
the study of 0v (3(3 decay. 
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